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ANNALS OF MATHEMATICS. 

Vol. III. February, 1887. No. i. 

A SPECIAL CASE OF THE LAPLACE CO-EFFICIENTS 3«. 

By Prof. Asaph Hall, Washington, D. C. 

(i.) When we consider two planets moving around the sun in the same plane, 
and take this plane for the plane of reference, the potential, or perturbative func- 
tion will be, in the notation of the Mecanique Celeste, 

m'r cos (z/' — v) m' 

r" v\^ + '"'^ — 2^^ cos (z^' — v)\ ' 

r, r', V, v' are the polar co-ordinates of the planets referred to the centre of the 
sun, and m' is the mass of the disturbing planet. If the eccentricities of the or- 
bits are small, the radii vectores will not differ much from the semi-major axes of 
the orbits. Let these be a and a' , and suppose a <^a' , a-=. a\a' ; also let v' — v 
■=■ d. Then we have to develop into a periodical series the function, 

I 

1/(1 — 2 a cos d ■+ 0?)' 

Omit the subscript s of Laplace, and assume 

== ^^(»> -I- d<-'> cos d + M^^ cos 2d + d^^> cos 3<? -f . . . 



1/(1 — 2« COS 0-\-(lf) 

+ <^"> COS td + . . . 

In the Mecanique Celeste the values of the 3'*' co-efficients are given in series 
which proceed according to powers of a. For the principal planets of our solar 
system this method is generally convenient, since the greatest value of a in this 
system is 0.723, which occurs in the theory of Venus and the Earth. But for 
values of a greater than this the series converge slowly, and its seems better to 
determine these co-efficients by a different process. In the second volume of his 
great work on elliptic integrals, Legendre shows how the ^<'' co-efficients may be 
reduced to these integrals, and their numerical values computed by means of his 
tables. In this case, however, the integrals can be reduced to forms so conve- 
nient that it seems worth while to proceed directly, and without reference to tables. 
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If we integrate the expression cos id cos JO dd, between :he limits o and n, 
we have, 



dd 



I COS id cosjd dd ^^\ I cos {i -\-j)d-\- cos {i —j) 
o o 

= ^r!!£ii±4^ + sin^WH] 

o 
Since i and/ are whole numbers, the integral becomes zero, except for i ^j, when 
its value is ^n. If i and / are both zero, the value of the integral is n. Denot- 
ing, therefore, the left-hand side of the equation for the series hyf{d), then mul- 

2 

tiplying the equation by - cos id dd, and integrating from o to tt, the terms are 
all zero except 

^''^ X - I cos id cos id dd. 
o 

We have, therefore, for the general value of the co-efficient, 

TT 

M'>=^f/{d) COS id dd. (i) 

o 

Thus for i = o, and ? = i, we have 

TT 

nJ t/(i — 2a cos d -\- a') ' 
o 

77 

,(j) .2 r COS d dd 

nJ |/(i — 2a cos d + a') ' 
o 

Let us consider three successive co-efficients of the series, 

TT 

^«- 1) = !//(<?) cos (i —i)d dd, 

O 

It 

b^'^ = ^f/(0) cos id dd, 

O 

b<' + ') = ^ //(<?) cos {i +i)ddd; 
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we can easily find a linear relation between these, co-efficients. Integrating by 
parts the expression 

cos id \/{l —2a cos d -\- a?), dd, 
we have 

I cos id 1/(1 — 2a cos 6 -\- a?) .dd 

sin id ,, a , 2\ « r sin ^ sin id dd 

= — ; — r/( I — 2a COS ^ -f- a^) — - . I -yj ^ 1 -ji\ • 

« "^ ^ ^ t J 1/(1 — 2a cos ^ + 0^) 

The first term is zero at the limits o and it, and since 

sin d sin id = ^ [cos (« — i) ^ — cos (? -f i) ^] , 

the result is 

f COS id 1/(1 — 2a cos d + oF)dd=— '^. M*-'^ + ^ . M' + '\ 
o 

Treating the same integral in a different manner, we have 
j cos id i/(i — 2a cos d + a") dd 

/ I ^\ r cos id dd r cos d cos id dd 

~ ^' + -"J |/(l - 2a cos d-]-a')~ ^V 1/(1 — 2a cos d + 0?) 

_, ^r cos id dd r cos (i + i)ddd 

^ "^ U 1/(1— 2a cos (? + «=') *'J ;/(i— 2acos^-i-a2) 

_ r cos (/ — i)d dd 
J \/{i — 2acosd + cF)' 

By comparison of the two values we have 

2« 2J V I y . 

or ^(i + i) = _^.i_±^.^«_£(zil.^(-i). (2) 

2? -f I a 2t -\- 1 ^ ' 

This equation shows that if we know the co-efficients ^<*~'> and ^<'', we can 
compute from their values the next co-efficient b^* + '\ Hence all the co-efficients 
may be computed from b'-"^ and ^''l We can find also from the definite integrals 
the derivatives of the l^'^ co-efficients with respect to a, which are required in the 
development of the disturbing forces, and these derivatives may be made to de- 
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pend on the ^**' co-efficients ; but these formulae are well known and appear to 
be convenient for computation. 

(2.) In order that we may reduce the expressions for 6^"^ and ^''' to forms 
suitable for numerical computation, we have to transform the integral containing 
j/(i — 2a cos d + (/) as a. divisor. Since cos d ^ 2 cos^ ^0 — i, we have 

|/(l — 2a cos d -\- a^)^ 1/(1 + 2a -\- 0? — 4a cos^ ^d) 
If, therefore, we put 



= (l+ 0)^(1 -^^^^cos^d). 



p - 4« „ _ 1 - t/(i - ^) 

""-{l+af «- i+i/(i_;P)' 

X ^ cos ■!■/?, dd =^ 77 ST ; 

^ ' i/(i — -^^ ) 

the value of d'-"^ becomes 

I 
m - 4 f d_^ 

-n{i-\-a)Jy^[ii-x^{i-Px^-]- 
o 

This elliptic integral may be transformed according to the method given by 
Lagrange.* For this purpose put 



^ = ''-^J[t^]- (3) 



Then ■^{i—Px')=i/[{i—x'){l- /^x')] , 

, dx X dx 

and 



V/[(I - x'') (I - /Px^)] -J/ (I - Bx') ■ 

Equation (3) gives f — J^x'''/ ^ x^ — x*, 

and differentiating, we have first, 

X dx dy 

(i — /^x'^y ~ i+J^y' — 2x^- 

Again, solving (3) lor x^, we find 

l+^f — 2X^ = -t/[l+2{^ — 2)f + k'y*2 , 

dy dy 



so that 



i+j^y' — 2x''~ ,/[i + 2 (/^ — 2) / + /^V] • 

*Oeuvres de Lagrange, Tome II. p. 253. 
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Assume the trinomial under the radical to be of the form (t + of) (i -|- /3/), 
and we shall have 

a = P — 2^ 2v/(i —k^ 

;9 = i^ — 2 — 2v/(l -/^); 

and the trinomial is resolved into the real factors, 

i + [i^-2 + 2v/(i-P)]y, 

I +[>^ — 2- 2V\l—P)']f. 

The quantities a and /3 are of the same sign since their product is k*^. The half 

sum of these quantities is .^ — 2, which is negative, and therefore a and y9 are 

negative. 

Let f = 2-P + 2y'{l-i^, 

^ = 2—I^ — 2y\\—ky, 

from which we have 

/=i + v/(i-n ?- 1-1/(1 -n 

Thus we have 

dx dy 



Now we have 



t/[(i - ;. ^) (I - /F^^)] ~ i/[(i -ff) (I - <ff)'\ 
2 2a 



and, if in equation (3) wc change our unit, or put 



our integral becomes 



/dx 
W\U-x')U 



I +« 



l/C( 



')(i-/^^^)] 



/ 



dy 



r[(i-/)(i-«!>'^)] 



This last integral is of the same form as the first one, and it is plain that we 
may continue such transformations as far as we please. Hence if we put as be- 
fore 

1/(1-^) _ 



+ 1/(1-^) 



- l/(i - <^ ) _ 



+ l/(i - ./) 

-j/(i - «,0 _ 
+ v/(i-«.') [I + i/(i 
etc., etc.; 



[I + y-(i - 


- ^)j^ ' 


(f 




[I + 1/(1 - 


-'^^)J'' 


«■' 





«I^)f 
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we shall have 

dx 

^-[(i-^^)(i-/F;.^)] 

^ I +« I + «i I +«2 I + «3 I +«.. r ^ 

2 ■ 2 '2 '2 ■ ■ ■ 2 V |/[(l — /)(l — a„y)]" 

From the values of a, a^, a^, . . . a„, we see that «„ can be diminished at 
pleasure, and, at last, a„y may be neglected. In this case we have 



Jwi-y)-*^- 



From equation (3) we see that y is zero when x ^o, and when x ^ \ ; and 
thatjy has a maximum for an intermediate value of x. Hence in the first step of 
the transformation we must double the integral for the limits, y ^oto y ^ \; 
and the same for each succeeding step ; so that I'n the final value we shall have 
the factor 2". The result will be 

^m ^ 4 I +« I +«1 I + <h I +«3 I + «n ^ 2» 

7t{i -\- a)' 2 ' 2 ■ 2 ■ 2 ••• 2 ••' 

or ^C) = 2(1+ a,) (i + o^) (I + a^) {i + a,) . . . (i + a„). (4) 

The value of b^^^ may be found by means of the same substitutions, which 
give 



m, ■ 4 r (2-y' —l).dx 

" — ;r(i +«)J vTd — -r2)(i— y^j 



■(I+«)J ^[(I_X^)(I_^;,^)]- 
o 

The preceding reductions give 

2x^ — I = /^y — p/[n- 2 (/P — 2) / + iy] ; 

and making the same change of the unit as before, we have 

2x^-x= \a{2f- I) + i« - i/[(i -f) (I - ^f)\. 
Putting, therefore, for dx\ j/[(i — •^^) (i — -^^^0]* its value found before, we have 

I II 

/(2;tr^ — \)dx \ -\- a a r {zf — \) dy \ -\- a r 

ySS}-x'){i-I^x^)\ - ""2~ • 2 J i/[(i-y)(i-«!;'^)] ^ J "^-^ 



o 



I 



+ 



\ -\- a a r dy 



J iTTCi 



2 • 2 -Ji/Ki-yjii-r//)]' 
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It is evident that this method of transformation may also be continued, and 
at last we shall have the forms 



/(2x^ — i)dx 



_ I 4 « I + «1 I + «2 

2 ■ 2 ■ 2 



2 ■ 2 ■ 2 ■ 2 



2 -J^/Ci-y) 



o o 

a ■ (I + a) (i + gQ (i + «,) . . . (i + a„) ^ 



-» + i 



.a,.(i +a)(i -l-ai) . . . (i + «„) 



+ < 



»+2 



a.ai.a,.(i + a)(i + aQ . 
-t- 2" + ^ 



(i + «„) 



i/y 



•J r/(i- 



/)• 



The first term approaches zero as its limit and may be neglected. The integrals 
of the second set vanish at the limits zero and unity. We have left only the last 
terms within the brackets, and since as before the value of the last integral is 
\Tt . 2", we have finally, 

^w = (i +«,) (1+ «,)(!+ a,) . . .(i+a„)|^a + -i + _0 + --^+ ...J , 



<5<" 



^^,<o,j^«^^. + «_^ + ^.A_M+. . .]. (5) 



Equations (4) and (5) can be computed easily, and then the remaining co- 
efficients may be found from equation (2), when a is not very small. 

We may notice here a relation between the moduli of the successive trans- 
formations which was pointed out by Lagrange, and which is similar to that em- 
ployed by Gauss in his celebrated memoir on the secular perturbations of a 
planet: Deierminatio Attractionis, etc. 1818. 

From the values of a, ai, a^, etc. we have 



_ ?! „ - 2? 

Pi Pi. 

and hence / is greater than q, pi greater than q^, etc. 



a = - , a, ^ — , a,^ — , etc. : 

/ A ' A 
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We have also /, = / + V{p'' - <f), q^=p- V{f - f) l 

p2=Pi + v\P' - q^\ ?2 = A - V{p' - qf) ; 

etc., etc. 

If we put p -\- q = m, pi -\- qi ^ nt^ ; 

p~q=n, pi — qi^fti; 
etc., etc. ; 

so that / = ^{m + «), q = ^{m — n) ; 

A = U^i + »i), q\ = i(»'i — «i) ; 

we shall have mi ^ m -\- n, «i := 2|/(;««) ; 

m2^t»i -{- «i, «2 = 2i/(wi«i) ; 
etc., etc. 

In the series of quantities, therefore, m, ttti, m^, . . . \ n, tii, n^, . . . , the 
corresponding terms are the arithmetical and geometrical means of twice the 
terms that precede them. 

(3.) For the principal planets of our solar system there does not appear to 
be much need of departing from the methods commonly given for computing the 
^^'' co-efficients ; but in the Saturnian system we find a case where these methods 
are not convenient. The two satellites of Saturn, Titan and Hyperion, the largest 
and the smallest of this system, move in the same plane, and from the values of 
their mean distances from the centre of the planet we have the ratio 

« ^ 0.8250863. 

Hence the usual series for computing the i'-'^ co-efficients will converge slowly, 
and the preceding method by definite integrals is much easier. I will compute 
the values of some of the.se co-efficients for this value of a. For this purpose put 

a = sin (p, «! = sin <pi, a^ = sin f 2, etc. 

• give 



The expressions for «!, a^, 0-3, 



— cos tp ^ , , 
— i- = tan^ ic, 

+ cos (p '■ 

— cos <pi 



-f cos (pi 
— cos ^2 



= tan" yi, 
= tan^ 1^2. 



and hence 



+ cos ^2 
I + «! = sec^ \<p, 
I + «2 = sec^ \<pi, 
I + «3 = sec== \(fi, 



etc.; 



etc. 
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The computation of ^*"' and ^*'^ by formulae (4) and (5) may be arranged as fol- 
lows : — 

a ^ sin f = 0.8250863, 

log sin <p = 9.9164993760. 



Calculation of b'-"''. 



<F =5S°35'5o'.'o7 

\f = 27 47 55-035 
^, = 16 8 16.349 






log log 

tan 9.7219831665 sec 0.0532568941 
sin 9.4439663330 
4 8.174 tan 9.1515781796 sec 0.0043209827 
= I 9 5.8228 sin 8.3031563592 

o 34 32.9114 tan 8.0021702287 sec 0.0000219319 

o o 20.8337 sin 6.0043404574 

o o 10.4168 tan 5.7033093 sec 0.0000000005 







0.0575998092 






2 log sec = 0.1151996184 






log 2 = 0.3010299957 






log tf^"^ = 0.4162296141 






^<"' = 2.607531805 




Calculation 


of6^'\ 


log a 


= 9.9164993760 


a = 0.8250863000 


logai 


= 9.4439663330 




logi 


= 9.6989700043 


2nd term = 0.1146662773 


2nd term 


= 9.0594357133 




log a, 


= 8.3031563592 


3rd term = o.OOi 1522908 


log* 


= 9.6989700043 




3rd term 


= 7.0615620768 


4th term = 0.0000000582 


log as 


= 6.0043404574 


0.9409049263 


logi 


=: 9.6989700043 


log = 9.9735457423 


4th term 


= 2.7648725385 


log i = 9-6989700043 


log a, 


= 1.4066186 


log b^"'' = 0.4162296141 


logi 


=: 9.6989700 


log <^<') = 0.0887453607 


5 th term 


= 3 870461 1 - 20 


/P^ = 1.226719759 



lO 
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Having found b'-"'' and b^^\ we can compute the remaining co-efficients by- 
formula (2). Thus 



(5<2); 



A(S) — 



= 1- 



I + O? 



I + a^ 



hm _ 1 . ^(0)_ 



bw _ I . bw^ 



' a ' 

etc., etc. 

The preceding calculations would be much shorter if made with seven-figure 
logarithms, and the accuracy would be sufficient for practical purposes ; but the 
work has been done with ten-figure logarithms in order to show better the ra- 
pidity of the convergence. The value of ^<"' computed by the common series, 
including twenty six terms and a'*, is 

^<'" = 2.60752853. 

The following table contains the values of the first six co efficients com- 
puted by means of the definite integrals and formula (2) ; and also the values of 
the same co-efficients computed from the formulae of Leverrier.* 

Values of (5j<'> ; a = 0.8250863. 



Co efficient. 


Value by preced-' 
ing method. 


Value by se:ies. 


^(0) 


2.6075318 


2.6072929 


^<'> 


1.2267 198 


1.2265599 


^(2) 


0.7967742 


0.7965316 


b<?) 


0.5624428 


0.5622774 


^a) 


0.4129396 


0.4126906 


6(5) 


O.3IO2719 


0.3100983 



The values computed by Leverrier's formulae include twelve terms, but they 
are too small, and it would be necessary to continue the series to double this 
number of terms if we wish to obtain values accurate in the sixth decimal place. 
But if series are to be used, they should be tran.sformed so as to proceed accord- 
ing to powers of i — a^ as was pointed out by Euler in 1750. 

(4.) The history of these co-efficients is interesting. Since it is necessary 



*Annales de P Observatoire de Paris, Tome II. Addition premiere, p. 5. 
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to consider the reciprocal of the distance between two bodies when their mutual 
attractions are to be determined, and the best practicable method of finding these 
attractions is generally the expansion of this reciprocal into a periodical series, 
mathematicians were immediately confronted with the problem of finding the co- 
efficients of the sines and cosines of this series. D'Alembert i^Jean le Rond) saw 
in 1743 that all the ^<'> co-efficients could be found as soon as the first two were 
known, a relation which is expressed in equation (2). Euler and Clairaut made 
many improvements in the actual numerical determination of these co-efficients. 
The elegant method of resolving the general expression of the reciprocal, 
(i — 2a cos ^ + «^)~', into two binomial factors by means of imaginary quanti- 
ties, which is employed by Laplace, is due to Lagrange. It is probable that the 
transformation of the elliptic integral which has been used above was suggested to 
Lagrange by the work of Landen. It has since been employed by several writers, 
and among others by Professor Woodhouse.* 

In looking at the great amount of labor that has been done on the impor- 
tant problem of developing the perturbative function, one cannot fail to notice the 
variety of the symbols that have been introduced. This variety has become so 
great as to cause confusion, and in order to be sure of his formulae it is now al- 
most necessary that every astronomer should have a development of his own. 
This condition seems to be unfortunate, since it causes a repetition of work that 
should be done once for all. 

U. S. N. Observatory, Nov. 5, 1886. 

Communicated by the Superintendent. 



ON THE FORM AND POSITION OF THE SEA-LEVEL AS DEPENDENT ON 
SUPERFICIAL MASSES SYMMETRICALLY DISPOSED WITH RE- 
SPECT TO A RADIUS OF THE EARTH'S SURFACE. 
By Mr. R. S. Woodward, Washington, D. C. 

[continued from vol. II, PAGE 1 3 1.] 

i8. We proceed now to the determination of the constants F„ of equation 
(6) and U^ of (56) and (57). It has already been stated that these constants are 
to be determined from the condition of equality in volumes contained by the dis- 
turbed and undisturbed surfaces, a condition whose analytical statem.ent is 

2~r^ i ri sin « da -=■ O. 
o 

* Philosophical Transaclioiis^ London, 1804. 



